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Abstract 

Nonmodal transient growth studies and estimation of optimal perturbations have been 
made for the compressible plane Couette flow with three-dimensional disturbances. The steady 
mean flow is characterized by a non-uniform shear-rate and a varying temperature across the 
wall-normal direction for an appropriate perfect gas model. The maximum amplification of 
perturbation energy over time, Gmax> is found to increase with increasing Reynolds number 
Re, but decreases with increasing Mach number M. More specifically, the optimal energy 
amplification Gopt (the supremum of Gmax over both the streamwise and spanwise wavenum- 
bers) is maximum in the incompressible limit and decreases monotonically as M increases. 
The corresponding optimal streamwise wavenumber, aopt, is non-zero at M = 0, increases 
with increasing M, reaching a maximum for some value of M and then decreases, eventually 
becoming zero at high Mach numbers. While the pure streamwise vortices are the optimal 
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patterns at high Mach numbers (in contrast to incompressible Couette flow), the modulated 
streamwise vortices are the optimal patterns for low-to-moderate values of the Mach number. 
Unlike in incompressible shear flows, the streamwise-independent modes in the present flow 
do not follow the scaling law G{t/ Re) ~ Re^, the reasons for which are shown to be tied 
to the dominance of some terms (related to density and temperature fluctuations) in the linear 
stability operator. Based on a detailed nonmodal energy analysis, we show that the transient 
energy growth occurs due to the transfer of energy from the mean flow to perturbations via an 
inviscid algebraic instability. The decrease of transient growth with increasing Mach number 
is also shown to be tied to the decrease in the energy transferred from the mean flow {£{) in 
the same limit. The sharp decay of the viscous eigenfunctions with increasing Mach number 
is responsible for the decrease of £i for the present mean flow. 



1 Introduction 

The linear stability (LS) analysis of compressible flows is of interest from the viewpoint of basic 
understanding on transition scenarios in such fluids [HI El El, and also because of its relevance in 
many high speed aerodynamic design problems. Since the growth of a small disturbance can lead 
to the flow breakdown, an understanding of different growth mechanisms becomes more relevant. 
The normal mode approach, where the initial value problem is reduced to an eigenvalue problem 
by considering the disturbance growth/decay in terms of an exponential time dependence, has 
been widely studied in the past. For more than a decade, however, the traditional modal stability 
analysis has been complemented by the analysis of transiently growing perturbations for various 
shear flow configurations |l4l|5l[6l|7l[8l. Such transient growth analyses have revealed that a flow 
can sustain large transient growth of energy in the parameter space that is stable according to the 
LS theory. It has been shown that the transient growth analysis can provide a possible reason for the 
experimentally observed critical Reynolds number (Recr) of various canonical flow configurations 
being less than that predicted by the LS theory (for a detailed review see the recent book of Schmid 
and Henningson ll9l). 
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It is well known that there are flow configurations (e.g. the plane Couette flow, pipe flow, etc.) 
that are stable according to the LS theory, but have been shown to have a finite Rccr in experiments. 
Therefore the nonlinear effects may destabilize these flows in the stable parameter region of the 
LS theory. However, for the nonlinearities in the governing equations to take over, the amplitudes 
of the disturbances must be large. Therefore, there may be a linear mechanism that causes an 
infinitesimally small disturbance already present in the flow to grow initially, so that the nonlinear- 
ities could act upon them. Such a mechanism for the transient growth exists due to the nonmodal 
growth caused by the non-orthogonality of the LS operator dH. In the stable region of the parame- 
ter space, even though each eigenstate decays during its evolution in space or time, a superposition 
of such eigenstates has potential for the transient inviscid growth before being stabilized by the 
viscosity at the rate of the least decaying mode. The experimental confirmation of this transition 
scenario comes from the observation that streaks precede the breakdown of laminar shear flows. 
Since the transient energy growth is maximum for streamwise vortices which subsequently give 
birth to streaks, it suggests that this transient growth mechanism provides a robust route to flow 
breakdown via disturbances that are stable according to LS theory. 

Even though the transient growth analysis has been extensively studied for incompressible 
flows E [lOl [6l El H m [m, similar studies on compressible flows are scarce lfT2l [T3l [T4l [Bll . 
For compressible boundary layers, Hanifi, Schmid and Henningson [fT2l have found the maximum 
transient energy growth {Gmax) increases with both Reynolds number {Re) and Mach number (M). 
They have used an energy norm such that the pressure-related energy transfer rate is equal to zero. 
They also cofirmed that the streamwise vortices are the optimal patterns as in incompressible flows. 
Farrell and loannou [14J have investigated the transient energy growth of 2D perturbations in the 
Couette flow (uniform shear flow) of a polytropic fluid with constant viscosity coefficients. They 
found that Gmax increases with Mach number. Recently, Tumin and Reshotko \15\ have developed 
a 'spatial' model of transient growth phenomena in compressible boundary layers. 

The present work deals with the 'temporal' stability analyses of the compressible plane Couette 
flow. The modal linear stability analysis of this flow has been reported by many investigators 
in the past [|T6l [TTl [TSl . In the incompressible limit this flow remains stable but instability is 
possible for a range of supersonic Mach numbers, the origin of which is tied to a family of acoustic 
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modes ifTTl [TSj. Duck, Erlebacher and Hussaini [ fTTjl have classified tliese acoustic modes into 
two distinct families depending on their phase speeds (see Fig. 2, the detailed description of 
this figure is provided in Sec. 3.1): the odd-modes (I, III, ...) that have phase speeds greater 
than unity in the limit of zero streamwise wavenumber (a), and the even-modes (II, IV, ...) that 
have phase speeds less than zero as a 0. (For the plane Couette flow, the non-dimensional 
velocities of the walls are bounded between and 1 since the top-wall velocity has been used 
for non-dimensionalization, see Fig. 1.) From an asymptotic analysis of the inviscid equations, 
they showed that the mode-I remains neutrally stable, but the mode-II can become unstable in the 
inviscid limit for large a. Hu and Zhong [flSl studied the same viscous plane Couette flow problem 
numerically, and reported new viscous instabilities. In particular, they showed that the mode-I 
can become unstable at finite Reynolds numbers due to the effects of viscosity but the ranges of 
Mach numbers and wavenumbers are very narrow. They further showed that the viscosity plays 
a dual role of stabilizing/destabilizing the mode-II instability. The inviscid mode-II instabilty of 
Duck et al. that occurs at large a is weakly stabilizied by viscosity, but the same mode-II is also 
unstable for low streamwise wavenumbers purely due to viscous effects [fTSl . Even though the 
higher-order even modes (IV,...) can also become unstable, the mode-II instablity represents the 
dominant instability for the compressible plane Couette flow. 

In this paper, the nonmodal transient energy growth analysis is reported for the compressible 
plane Couette flow of an appropriate perfect gas model for three-dimensional disturbances. The 
governing equations and the mean flow are detailed in Sec. 2. The linear stability problem is 
formulated in Sec. 3, and the related numerical method and its validation are detailed in Sec. 3.1. 
For the tranisent growth analysis, the disturbance size is measured in terms of the Mack energy- 
norm BHIHl. The results on the energy growth, the optimal perturbations and their structure, and the 
scalings of maximum transient growth are presented in Sec. 4. The constituent energies transferred 
from the mean flow, the viscous dissipation, the thermal diffusion and the shear-work during the 
transient growth are investigated in Sec. 5 for the initial perturbation configuration that would reach 
the optimal configuration at a later time. The inviscid limit of stability equations are analysed in 
Sec. 5.1. The conclusions are provided in Sec. 6. 
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2 Governing Equations 



We consider the plane Couette flow of a perfect gas of density p* and temperature T*, driven by 
the relative motion of two parallel walls that are separated by a distance h*. The top wall moves 
with a velocity and the lower wall is stationary, with the top-wall temperature being held fixed 
at T^. Note that the quantities with a superscript * are dimensional, and the subscript 1 refers to 
the quantities at the top wall. The velocity field is denoted by u*(x*, y*, z*,t*) = {u*, v*, w*)"^, 
with u*, V* and w* being the velocity components in the streamwise (x*), wall-normal (y*) and 
spanwise (z*) directions, respectively. For non-dimensionalization, we use the separation between 
the two walls h* as the length scale, the top wall velocity, f/j", and temperature, T^, as the velocity 
and temperature scale, respectively, and the inverse of the overall shear-rate, U*/h*, as the time 
scale. The governing equations are the three momentum equations, continuity, energy and state 
equations which are omitted for the sake of brevity. The flow is described by the Reynolds number 
Re, the Prandtl number a and the Mach number M, defined via 

Re = U*plh*/i2l, a = ^i*Cp/K*, and M = U*/^^RT^. (1) 

Here 7 = c*/c* is the ratio of specific heats, K* the thermal conductivity, and R = c* — cl the 
universal gas constant; the Prandtl number is assumed to be a constant, a = 0.72 and 7 = 1.4. For 
the shear viscosity p, we use the Sutherland formula 

^= ^'(T + c^^ ' with C = 0.5. (2) 

Following Stokes' assumption, we set the bulk viscosity to zero (i.e., C = 0) such that 

2 2 

To make a direct comparison with earlier linear stability results of the compressible plane Couette 
flow ETIIIHI, we have used the same constitutive model for the transport coefficients of a perfect 
gas. 
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2.1 Mean flow 

The mean flow is that of the steady, fully developed plane Couette flow for which the mean fields 
are given by 

[u,v,w] = [Uo{y),0,0], p = Po{y), T = Toiy). (3) 

Hereafter the subscript is used to refer to the mean flow quantities. The continuity and the 
2-momentum equations are identically satisfied for this mean flow. From the y-momentum equa- 
tion, it is straightforward to verify that the mean pressure, po, is constant, chosen to be pq = 1. 
Hence, from the equation of state, the mean density is related to the mean temperature via po{y) = 
l/To(y). The remaining x-momentum and energy equations for the mean flow are 

which have to be solved, satisfying the following boundary conditions: 

f/o(0)=0, t/o(l) = l, To(0)=T^, ro(l) = l, (6) 

where T^, is the non-dimensional temperature of the lower wall. The energy equation can be solved 
exactly to yield the temparature in terms of the flow velocity [17] : 



To 

where 



r + {l-r)Uo-[l-^]U'^ 



(7) 



= 1 + [(7 - l)aM^]/2 (8) 

is the recovery temperature, and r = T^/Ty. is the temperature ratio. The streamwise velocity, 
Uo{y), is then calculated numerically by solving the following equation: 

^ = (9) 

dy piUo)' 

using a fourth-order Runge-Kutta method. The shear stress, r, is a constant that must be deter- 
mined iteratively such that Uo{y) satisfies its boundary values: 

f/o(0) = and Uo{l) = 1. (10) 



With r = 1, which corresponds to an adiabatic lower wall, the calculated profiles of the mean 
velocity (Uq), the temperature (Tq), the shear rate (Uoy) and the viscosity (/io) are shown in Fig.[T] 
for two representative values of the Mach number, M = 2 and 5. (Note that the mean flow profiles 
do not depend on Reynolds number.) Clearly, the shear-rate, the temperature and the viscosity are 
non-uniform along the wall-normal direction; the deviation from the corresponding uniform shear 
flow increases with increasing M. We shall return back to discuss the effects of non-uniform mean 
flow on certain aspects of the transient growth results in Sec. 5.1. 



3 Linear Perturbation System 

We impose three-dimensional (3D) perturbations on the mean flow described in Sec. 12.11 The flow 
quantities with small- amplitude disturbances are considered as q = qo + q, where qo represents 
any mean flow quantity, and q is its perturbation. For qo ^ q, the governing equations are 
linearised with respect to q. Then the solutions of the linearized equations are assumed in the form 
of normal modes, i.e., 

ci{x,y,z,t) = ci{y,t) exp [i{ax + /3z)], (11) 
q(y,^) = q'(2/) exp (-zc^t), (12) 

where a and P are the streamwise and spanwise wavenumbers, respectively, and to = ujr+iuji is the 
complex frequency. The imaginary part of the complex frequency, tUj, represents the growth/decay 
rate of perturbations, and its real part, uj^, is related to the phase speed of the perturbation via 
Cr = ijJrjoL. The flow is said to be asymptotically stable if cuj < 0, unstable if cuj > and neutrally 
stable if uJi = 0. Later for the transient growth analysis (Sec. 4), a summation over all tu's with 
appropriate coefficients will be implied on the right-hand-side of equation (fT2l) . 

The perturbation pressure is removed from the linearized equations by using the perturbed 
equation of state 

p = np + f/To. (13) 
Similarly, the perturbation viscosity is written in terms of perturbation temperature: 

ft = firf, 

7 




Figure 1: Profiles of the mean velocity Uq, temperature Tq, shear-rate Uoy and viscosity //q for 
M = 2 ( ) and M = 5 ( ). 
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where fir = (d/u/dT)o is evaluated at mean flow conditions. FinaUy, the reduced ordinary differ- 
ential eigensystem in five unknowns can be written as 

£q' = cjlq', (14) 

where q' = {u' ,v' ,w' , p' ,T'}^ is the eigenf unction and I the identity matrix. The elements of the 
linear operator C are given in the appendix. The boundary conditions are: 

dT' 

u'{0) = u'{l) = 0; v'{0) = v'{l) = 0; w'{0) = w'{l) = 0; T'(l) = —(0) = 0. (15) 

Note that the temperature boundary condition corresponds to an isothermal upper wall and an 
adiabatic lower wall as in the work of Duck et al. ifTTl and Hu and Zhong [fTSl . 

3.1 Numerical method and code validation 

The differential eigenvalue problem (fT4l)-(fT5l) is discretized by using the Chebyshev spectral method |fT9l . 
The transformation ^ = {2y — 1) is used to map the physical domain y G [0, 1] to the Chebyshev 
domain ^ G [—1, 1]. The equations are collocated at the Gauss-Lobatto collocation points: 

e, = cos(^), j = 0,l,2,...,iV, (16) 

which are the extrema of the A^-order Chebyshev polynomial Tjy. There are two sub-methods of 
the Chebyshev spectral method to discretize this eigenvalue problem (fT4l) . In one method, any 
function 0, say, is approximated by using a Chebyshev polynomial of degree A^: 

</)(0 = i:anT„(0, (17) 

n=0 

with the Chebyshev coefficients a„ being treated as 'unknowns'. The derivatives of 0(^) are also 
found in terms of these polynomials through a recurrence relation that they satisfy [9J. In the 
other method, the function is expressed in terms of its values at the collocation points using an 
interpolant that uses Chebyshev polynomials, and the derivatives are found from an interpolation 
formula (see Malik [|19i for details). This latter method is convenient since it results in a simple 
eigensystem, whereas the former results in a generalized eigensystem. The x, y, 2;-momentum 
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equations and the energy equation are replaced by the boundary conditions for u' , v' , w' and T' at 
walls, respectively. The unused y-momentum equation is used to replace the continuity equation 
at the walls (also termed as artificial boundary condition [19]). Since the boundary conditions are 
independent of the eigenvalue cu, the right hand side of the eqn (fT4l) would be singular, as the rows 
corresponding to boundary points would be zero. This singularity is removed by the row/column 
operations, resulting in a reduced system. The QR-algorithm of the Matlab software is used to 
solve this eigenvalue system. 

To validate the stability code, we have compared our numerical results on asymptotic stability 
with the results of Hu and Zhong [fTSll . For example, for the parameter values of -Re = 2 x 10^, 
M = 2, a = 0.1 and P = with N = 100 collocation points, our data on both growth rate and 
frequency match those of Hu and Zhong upto the fourth decimal place. For this parameter set, the 
distribution of eigenvalues, c = to /a, is shown in Fig. [2tct), and the enlargement of Fig. WLcl) is 
shown in Fig.|2&) that portrays the well-known 'Y'-branch of the spectra that belongs to viscous 
modes. Following the classification of Duck et al. [[TtI . the odd-family (I, III, ...) and even-family 
(II, IV, ...) of inviscid/acoustic modes are shown in Fig.^a); the mode-I and mode-II are indicated 
in Fig.[2t6). Note that the 'Y' -spectrum for compressible Couette flow becomes difficult to resolve 
numerically if Re and a are large ifTTlfTSl . 

Now we present a few representative stability diagrams for compressible plane Couette flow, 
delineating the zones of stability and instability. Figure IS^a) shows the contours of the growth rate 
of the least decaying mode, cuidi = max(c<jj), in the (Re, a)-plane for two-dimensional disturbances 
(P = 0) at a Mach number of M = 5. (This plot matches well with Fig. 20a of Hu and Zhong ifTSl .) 
For 3D disturbances. Fig. ^b) shows the contours of coidi in the (a,/9)-plane for M = 5 and 
i?e = 2 X 10^. In each plot, the flow is unstable inside the neutral contour, denoted by '0', and 
stable outside. The mode-II (see Fig. 2b) is responsible for the observed instabilities; for other 
details on linear stability results the reader is referred to Hu and Zhong [fTSl . 

For the transient growth analysis, we will focus on the control parameter space (such as in 
Fig. 3) where the flow is stable according to the modal linear stability analysis, and investigate the 
potential of this 'stable' flow to give rise to transient energy growth. 
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Figure 2: Spectra for M = 2 and i^e = 2 x 10^: a = 0.1 = = 100. Panel h is the zoom of 
panel a. 
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Figure 3: (a) Contours of the growth rate of the least-decaying mode, uidi, in the (Re, Q;)-plane 
with M = 5 for two-dimensional (/3 = 0) disturbances, (b) Contours of coidi in the (a, /3)-plane 
with Re = 2x 10^ and M = 5. 
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4 Transient Energy Growth 



To study the evolution of a general non-eigenstate, we relax the assumption of exponential time 
dependence. The linear perturbation equations can be written as [HI 

where q(?/, t; a, /3) is the inverse Fourier transform of q(x, y, z, t), as given by (fTTI) . The state 
function is expanded in the basis of its eigenfunctions as 

q=E'^4^)qfc> (19) 

k 

where the index k runs through a selected portion of the least decaying eigenmodes in the complex 
eigenvalue plane (refer to Fig. 2). Using (fT9]) . the equation (fTSl) can be diagonalized as 

— = -zQk, with K{t) = e-''^K{0), (20) 

(J.6 

where k = {nk}'^ and = diag{u;fc}. 

To compute the nonmodal energy growth, we need an expression for the perturbation energy. 
Under a suitable definition of the inner product, the square of the norm ||q(t)|| can be made to 
represent a measure of this energy density, £{a,P,t). The energy density is defined as that of 
Mack Em 

S{a, P, t) = qt(e, t)Mq,{i, t)de, (21) 

where the weight matrix, 

M = diag{po,Po,Po,To/po7M^po/7(7 - 1)ToM2}, (22) 

is diagonal. In this definition of S{a, P, t), the spatial average of the rate of pressure-related work 
(i.e. the compression work) is zero. As elaborated by Mack [[3l, the contribution of this com- 
pression work to the total perturbation energy should vanish due to the conservative nature of 
the compression work. Note that the Mach energy norm (|2TI) has also been used in the 'spatial' 
transient growth analyses of compressible flows IfTSll . 
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Let G{t,a, /3; Re, M) be the amplification of the initial energy density maximized over all 
possible initial conditions, i.e., 

Git, a, /3; Re, M)^ Git) = max|f^lM. (23) 

q(o) t{a,fJ,{)) 

For the numerical evaluation of Git), we refer the readers to Refs. [[6ll9l[T2|. 

A typical variation of the energy amplification factor. Git, a, (3), with time is shown in Fig.llfa), 
for the parameter values of a = 0.1 = i?e = 2 x 10^ and M = 2. There is a 200-fold increase of 
the initial energy density at t ^ 300. Since the flow is stable for this parameter combination, G(t) 
decays in the asymptotic limit (t oo). To compute Git) we do not need to consider all modes 
since we have found that the modes with large phase speeds (see Fig. 2) do not contribute to the 
energy growth. This observation is similar to the earlier findings for compressible boundary lay- 
ers [[T2|. Therefore, to reduce the computational effort, the index k in (fT9l) is chosen corresponding 
to the modes whose phase speeds are within the range — 1 < uj^/a < 2 (i.e., comparable to the 
extremes of the mean flow velocity which varies between and 1), and the decay rate is less than 
0.5 (i.e., uji > —0.5). With this choice of modes, the maximum possible energy amplification over 
time, 

Gmax(a,/5;^e,M) = maxG'(t,a,/3;i?e,M) at t = t^ax, (24) 

saturates to some constant value with the number of modes retained K, as shown in Fig.|56). It is 
observed that out of a total of 500 modes, only 70 modes are sufficient to compute Gmax accurately, 
and the contribution of the remaining modes to Gmax is negligibly small. 

For M = 3 and /3 = 3, the contours of Gmax in the (Re, Q;)-plane are shown in Fig. I^a). For 
other values of M and P, the Gmax -contours look similar. It is observed that Gmax varies non- 
monotonically with streamwise wavenumber a for fixed Re, and is maximum for non-zero values 
of a; Gmax increases with Reynolds number for any value of a. The latter observation mirrors 
similar findings on transient growth in incompressible fluids flU. 

Moving onto the effect of Mach number on the transient energy growth, we show, in Fig. |6ta) 
and|6t6), the contours of Gmax in the (a,/3)-plane for M = 2 and 5, respectively; the Reynolds 
number is set to Re = 10^. It is observed that the magnitude of G^ax decreases sharply as we 
increase the Mach number from 2 to 5 for any combination of a and (3. The supremum of Gmax 
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Figure 4: (a) Variation of the energy amplification factor, G{t), with time, (b) Dependence of the 
maximum energy growth, G^ax^ in panel a with the number of modes, K, considered (N — 100). 
Parameter values are i2e = 2 x 10^, M = 2 and a = 0.1 = jS. 
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Figure 5: Contours of maximum transient growth, Gmax, in the (i?e, Q;)-plane for M = 3 and 
/? = 3. 
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Figure 6: Contours of Gmax in the {a, /3)-plane for Re — 10^ at various Mach numbers: (a) M = 2; 

(b) M = 5. 
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over all possible combinations of wavenubers (a, P) is called the optimal energy growth, Gopt^ 

Gopt(i?e,M) = snpG^axia,/3;Re,M), (25) 

which corresponds to (topt, «opt, Popt)- In Fig- 13^, b), this global maximum is seen to occur at 
a ~ and /3 ~ 3. Also, Gopt decreases by an order of magnitude with increasing Mach number 
as seen in Figs. [6ta) and 06). (An explanation for this behaviour is provided in the next section, 
along with a detailed energy analysis in terms of different constituent energies.) It is worth pointing 
out that our result on the decrease of Gmax with M is in contrast to that for boundary layers [fT2]| . 

The variations of the optimal growth Gopt, the corresponding optimal time topt, the optimal 
streamwise, Oopt^ and transverse, Popt, wavenumbers with the Mach number are shown in Fig.|3a- 
d); as in Fig. 6, Re = 10^ for these plots. It is observed that while both Gopt and topt decrease 
monotonically with increasing M, both the optimal wavenumbers (ctopt and /3opt) vary nonmono- 
tonically in the same limit. (Note that the wiggles in panels c and d are due to the searching of 
Gopt in the {a, /3)-plane with small but finite grid-spacing- these computations are very time con- 
suming.) Figure Ulc) shows the interesting result that the optimal streamwise wavenumber Oopt 
can become zero at large M. This is in contrast to the incompressible Couette flow for which Oopt 
stays always slightly away from the null value [|91. 

4.1 Structure of optimal disturbances 

The initial disturbance pattern for given values of a and P that reaches the maximum possible 
transient growth Gmax at tmax along with the corresponding disturbance pattern at time tmax can 
be determined through the singular value decomposition [9|. In the following, such disturbances 
at t = and tmax are referred to as the optimal disturbances. 

Before presenting the structure of optimal velocity patterns, we show the variation ofS{t)/£(0) 
with time, denoted by the dash line, in Fig. [8] for the particular initial perturbation that reaches 
maximum energy, Gmax, at a time tmax- The parameter values are set to M = 2, Re = 2 x 10^, 
a = and P = 3. For comparison, we have also superimposed the variation of G{t) with time, 
denoted by the solid line. As expected, the envelope of £{t)/£{0) is below that of G(t) since the 
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Figure 7: Variations of (a) the optimal energy growth, Gopt, (b) optimal time, topt^ (c) optimal 
streamwise wavenumber, Q!opt> and {d) optimal span wise wavenumber, /3opt, with the Mach number 
M for Re = 10^ 
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Figure 8: Variations of G{i) (solid line) and £{t)/£{Q) (dash line) with time for a = 0, /3 = 3, 
M = 2 and i?e = 2 X 10^ 
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former corresponds to a particular initial condition but the latter to all possible combinations of 
initial conditions that yields maximum amplification at any time. 

For the maximum energy growth as in Fig. [8l the optimal velocity patterns at t = and tmax 
in the (y, z) -plane are shown in Fig.|9ta) and Fig. [9^6), respectively. (Note that these patterns are 
invariant along the streamwise direction since a = 0.) Both patterns resemble streamwise vortices 
which are typical of all shear flows. It is noteworthy that the optimal pattern in Fig. (9^6) has two 
counter-rotating streamwise vortices along the wall-normal direction. The inviscid rapid growth 
of these optimal patterns for null or very low values of the streamwise wavenumber (a) is due to 
the lift-up mechanism ||20| during the motion of the fluid particles in the wall-normal direction. 
Another typical feature of shear flows is the spanwise alternating streamwise velocities as shown 
inFig.Htc). 

Figures \Wia) and [T0l 6) show representative two-dimensional {(3 = 0) optimal velocity pat- 
terns in the (x, y) plane. The parameter values are set to « = 1, M = 2 and Re = 2 x 10^. It is 
observed that the pattern at t = has structures that locally oppose the mean shear. This pattern 
subsequently evolves into two 'cross-flow' vortices at t = tmax- This final configuration is the 
outcome of the well-known Orr-mechanism ^ that leads to transient energy growth due to the 
tilting of the initial perturbations into the direction of the mean shear as time evolves. 

Returning to the variation of the optimal wavenumber, aopt, with the Mach number in Fig.|7tc), 
we note that aopt is very small but finite at low-to-moderate values of M and zero at some large 
value of M. This implies that the pure streamwise vortices are the optimal patterns at high Mach 
numbers, but their modulated cousins are the optimal patterns for low-to-moderate values of the 
Mach number. Our results on optimal patterns at high Mach numbers (M > 5.5) should be con- 
trasted with that for incompressible Couette flow for which the oblique modes constitute optimal 
patterns. Typical optimal patterns of such modulated streamwise vortices are shown in the {x, y)- 
and {z, 2/)-plane in Fig. [TT] for the optimal energy growth (Gopt) at M = 5, with other parame- 
ters as in Fig. 7. We conclude that the structural features of the optimal disturbance patterns in 
compressible Couette flow look similar to those of incompressible shear flows. 
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Figure 9: Optimal patterns of perturbation velocity at (a) t = and (6) t = tmax in the (z, y)-p\ane 
for a — 0, P — 3, M — 2 and i?e 2 x 10^. (c) Contours of optimal streamwise perturbation 
velocity att — t^ax for the same parameter set. 
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Figure 10: Optimal velocity patterns at (a) t = and (b) t = tmax in the (x, 7/)-plane for two- 
dimensional disturbances {(3 — 0) with a — 1, M — 2 and Re — 2 x 10^. 
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Figure 11: Optimal velocity patterns at (a) t = and (b) t = t„iax in the (x,y)-p\ane. at z = 0; 
optimal velocity patterns at (c) t = Q and (d) t — t^ax in the (z, y)-p\ane at x = 0. Parameter 
values correspond to the optimal energy growth Gopt in Fig. 1: M — 5, a — 0.02, /3 = 3.5 and 
Re = 10^ 
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4.2 Scalings of G,nax and tmax 

For incompressible channel flows with streamwise-independent modes (i.e., a = 0), Gustavs- 
son [ITOl has shown that Gmax varies quadratically with the Reynolds number Re, and t^ax varies 
linearly with Re. More specifically, when the energy growth curve, such as the one in Fig. 4(a), is 
plotted in terms of G{t) / Re^ and t/Re, the renormalized growth curves for different Re collapse 
onto a single 'universal' curve. Following a similar analysis, Hanifi and Henningson lfT3l found 
that this scaling law also holds for compressible boundary layers. However, for the present Cou- 
ette flow of non-uniform shear and non-isothermal fluid, this scaling law does not hold as shown 
in Fig. [T2] which displays plots of G{t) / Re^ verses t/Re for a wide range of Re at M = 2. Similar 
trends persist at other values of Mach number (not shown for brevity). 

For the present flow configuration, it can be verified that the Mack transformation [|3l [131 

{u', v', w', p' , T', t} {Re u, V, w, Re p, Re T, Re t} (26) 

(hence p' — > Re p, from the equation of state), can make the streamwise-independent linearized 
stability equation (fTSi) independent of the Reynolds number Re, except for the terms associated 
with density and temperature fluctuations in the y- and 2;-momentum equations (i.e., ^25, ^34 
and £35 as detailed in the Appendix). Only on neglecting these terms, one can show that Gmax ~ 
Re^ as detailed by Hanifi and Henningson [fT3l who further assumed that p' p/ Re to neglect 
£24, £25, £34 and £35 for large Re. The above pressure-related terms may not be negligible for all 
mean flows. This can be ascertained for the present flow configuration if we recompute the energy 
growth by setting 

£24 = £25 = ^34 = £35 = (27) 

in the linear operator £. Indeed the rescaled growth curves for different Re now collapse onto a 
single curve as shown in Figs. \T3t a-b), for the same parameter values as in Figs. \TT[a-b). (Note that 
in these plots the energy growth does not decay since the above procedure (|27l) introduces some 
'artificial' neutral modes in £, resulting in an asymptotic value for G{t) at large times.) Therefore, 
we conclude that the non-negligible values of £24, £25, £34 and £35 for the plane Couette flow are 
responsible for the invalidity of the scaling law Gmax ~ Re^- 
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Figure 12: Verification of the quadratic scaling law for G{t) for streamwise-independent (a 
disturbances, (a) = 0.1 and M = 2; (b) /? = 1.0 and M = 2. 
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Figure 13: Verification of the quadratic scaling law for G{t) for streamwise-independent {a — 0) 
disturbances without pressure-related terms in C (see text for details), (a) /3 = 0.1 and M = 2; (b) 

(5 = 1.0 and M = 2. 
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We note in passing that Hanifi and Henningson's [fT3ll assumption of p' p/ Re leads to a 
change in the equation of state by introducing a factor of Re into it. The form-invariance of the 
equation of state is essential since it has been used to remove p' from the linear perturbation system. 
In general, the (streamwise-independent) linear operator C of compressible flows cannot be made 
free from the apearance of Re (via the Mack transformation), and hence the scaling, Gmax ~ -Re^, 
would not hold for all mean flows. For special cases, the terms £24, ^25, and £35 might be 
negligible (e.g. in compressible boundary layers [13J at high Re) and hence the scaling law would 
hold there. 

5 Nonmodal Energy Budget 

We now consider the evolution equation for the perturbation energy £(a, P,t). Multiplying (fTSi) by 
q^A^ and adding the complex conjugate of the resulting equation to itself, we obtain the evolution 
equation for £{a, P, t) as 



with c.c. being the complex-conjugate term. In the above equation, the total perturbation energy 




(28) 



has been decomposed into several constituent energies, Sf 




To7(7 - l)M- 



-T^t> dy + c.c. 



(29) 




+ (/xo + Xo)D)v + a/3{fio + Xo)u^ui 
-{UoyyPT + UoyTQyfiTT)u^T - UoyfiTU^DT 
-iav\Xoy + {fio + Xo)D)u 
+IJo{a^ + f3'^)v^v - v\{Xoy + Poy)D 
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+ (Ao + fio)D^ + fioyD + fioD'^)v 
— i/5(Ao + ^q)v'^ Dw — iaUoyfiTv'^T 
—i^XoyV^'w — i(3^QyW^v 
+a/?(/io + \q)w'^u — i/?(Ao + fiQjw'^ Dv 
+ (/io(a' + + /5'(Ao + ^io)Ww 
—IXqw^D^w — ^QyW^ Dw dy + c.c. 



(30) 



a/?e(7 - 1)M2 



(31) 
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-Re 7o 



Po 



>u 



+ 2iafioUoyT^v + Uly^rT'^T dy + c.c. 



(32) 



Here, ^1 represents the energy transfer from the mean flow, £2 is the viscous dissipation, £3 the 
thermal diffusion and f'4 the shear-work, respectively. Note that there is an additional term, £q. 



£, 



-%a 



P07M2 



p^p + 



Po 



ro7(7 - 1)M' 



0, 



d^ + c.c. 

(33) 



representing the convective transfer of perturbation energy (by the mean flow), which is identically 
zero. The decomposition of the total perturbation energy into different constituents, as in (I29H32I) . 
is useful to analyse the role of each constituent energy on the growth/decay of total energy. The 
evolution equation (l28l) provides an energy-budget for the total perturbation energy, and helps 
to quantify the contributions of different kinds of perturbation energies, leading to the transient 
growth. 

To analyse the nonmodal energy budget, we choose a set of values for a, /5, M and Re that 
leads to transient growth. The constituent energies, £j, are then each evaluated using a quadrature 
formula with an initial perturbation configuration that would reach the optimal energy. The initial 
values of the constituent energies are chosen as {£'(0)} = {1,0, 0, 0, 0} so that the total intitial 
energy is equal to the normalized value. 
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Figure 14: Energy budget with time for M = 2 and Re — 2 x 10^; solid line, total energy; 
dash line, transferred from mean flow; dot-dash line, lost viscous dissipation; marked line, lost by 
thermal diffusion, (a) a = 0.1, /3 = 0.1. (b) a = and /9 = 3 
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Figure [Eta) shows the energy budget for three-dimensional disturbances with a = 0.1 = /?; 
other parameters are as in Fig. 4(a). For this case, the energy gain by the shear-work term, f'4, is 
found to be negligible and hence not shown in this plot. It is observed that the energy transferred 
from the mean flow. Si, (denoted by the dash line) increases with time, reaching a large asymptotic 
value beyond the optimal time (= tmax)\ this is the major cause for the transient growth. At some 
later time {t > tmax), this transient growth is nullified by the energy taken away by the thermal 
diffusion, ^3, (denoted by the line with solid symbols) and the viscous dissipation, £2 (dot-dash 
line). For this parameter set, the energy loss due to the viscous dissipation dominates over that due 
to the thermal diffusion at large times. 

Figure UA^b) shows the same energy budget for a representative case of streamwise independent 
{a = 0) disturbances, with parameter values as in Fig. 8. As in the previous case (Fig. [T4b). the 
energy transferred from the mean flow is responsible for the transient growth; however, unlike in 
Fig. [T4T a). £1 reaches a peak and then decreases to attain an asymptotic value. The variations of 
the other constituent energies are similar to those in Fig. [TJa). 

In the previous section, we found that the optimal transient energy growth Gopt decreases with 
increasing Mach number. To explain this behaviour in terms of constituent energies, we show 
the energy budget versus Mach number in Fig. [15] for the optimal values of {a, P, t) (as in Fig. 
[7]). Note that the total energy curve shown in Fig. \T5\a) (denoted by the solid line) matches with 
the one in Fig. I3a). It is observed that the energy transferred from the mean flow (£1) and the 
viscous dissipation loss {£2) are the most significant compared to the other components {£3 and £4, 
see Fig. [TSb ) of the perturbation energy. Since £1 decreases with increasing M, so does the total 
energy at the optimal time. In general, the decrease in the transient energy growth with increasing 
M is primarily due to the decrease of the energy transfer from the mean flow to perturbations in 
the same limit. 

For the present flow configuration, it is still not clear why the transferred energy from the mean 
flow, £1, decreases with increasing Mach number; an inviscid energy analysis provides an answer 
for this as well as for the transient growth mechanism. 
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Figure 15: Energy budget with Mach number at (ctopt, /9opt, ^opt) for Re — 10^. (a) solid line, total 
energy; dash line, transferred from mean flow; dot-dash line, lost viscous dissipation, (b) solid 
line, energy lost by thermal diffusion; dash line, energy gained by the shear-work term. 
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5.1 Inviscid limit: Mechanism for transient growth 

For streamwise independent disturbances (a = 0), Hanifi and Henningson lfT3l have shown that 
the inviscid linear stability equations of compressible fluids have the following solution 

qivs(2/,t) = iqivs(2/,0), (34) 

where 

A = l + A and qivs(l/, t) = {mivs, t'ivs, Pivs, ^Ivs}'^, (35) 

and the subscript, 'ivs' stands for 'inviscid'. The nonzero elements of the matrix, A, are A12 = 
—Uoyt, A32 = —poyt and ^442 = —Toyt. In fact, this inviscid solution for velocity perturbations is 
the same as that of EUingsen and Palm's incompressible solution [|2T|. The inviscid perturbation 
energy, Gi^s{t), maximized over q-„s{y, 0), in the same definition of the energy norm, is given by 

Givslt) = max({Afc}), (36) 

where {A^j's are the eigenvalues of the differential equation 

A^CAqi^siy, 0) = XCqi^y, 0), (37) 

with C being the associated linear differential operator 

C = diag{po,Po(l -r'd/dy),roV7M2,p2/^(^ _ i)^^}. (38) 

The above equation is solved using the same spectral method described in Sec. 3.1 with the fol- 
lowing boundary conditions: 

i)i,s(0,0) = 5i,s(l,0) =0. (39) 

Figure [16] shows the variation of G'ivs(i) with time for P = 1 and M = 2. For comparison, 
the viscous G'(t) -curves at various Re are also displayed. In terms of different energy components, 
the inviscid energy growth rate of the optimal perturbation configuration is given by (|29|) . since 
the transfer rates of the other energy components, £2, £3 and £4 are zero in the inviscid limit. 
Therefore the reason for the inviscid growth is due to the energy transferred from the mean flow to 
perturbations, i.e. £1. This energy transfer occurs due to an algebraic instability [l20l|2T|, wherein 
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the streamwise perturbation velocity, density and temperature grow algebraically (linearly) with 
time, if the initial normal perturbation velocity is non-zero (i.e. Vi^siO) 7^ 0). A noteworthy point 
in Fig. [T6]that the inviscid growth-curve, G'ivs(i), coincides with the viscous curves only /or a very 
short time- we shall return to explain this point later. Note that this observation is in contrast to 
boundary layers (see Figs. 1 and 2 in ref. [fT3l ) for which the viscous growth curves (at different 
Re) closely follow the inviscid growth-curve till they achieve their maxima Gmax- 

The inset in Fig. [16] shows the curves of Givs{t) (= £i) for three different Mach numbers. 
Clearly, at any time, £i increases with increasing M for the inviscid case. This is in contrast to the 
viscous case where we have seen a steady decrease of £i with increasing M (dashed line in Fig. 
15a). To understand this difference, let us rewrite the energy transfer rate (|29l) as 

^i(t) = - q^MqoyV dy + c.c, (40) 

which is identical for both the inviscid and viscous cases. Now combining (|40|) with (l34l) . we 
obtain 

£Ut) = -2 £ qL(l/, 0)A^CqoyvUy, 0) dy (41) 

In this equation the M-dependence is present only via A'^Cqoy Note that the inviscid temperature 
eigenfunction, 

fi^sit) = TiysiO) - Viysi0)T0yt, (42) 

grows quadratically with M since T^y ~ M^. It can be verified that the dominant contribution 
to £ivs{t) in (|4T)) at large M comes from the energy associated with temperature fluctuations, and 
hence, for any given initial condition, the norm of £nsit) would increase with increasing M. We 
can conclude that the increase of the inviscid energy growth with increasing M is primarily due to 
the increased energy transfer from the mean flow to temperature fluctuations. 

Though the above contribution is also present for the viscous case, the decay of the viscous 
eigenfunction q by viscosity becomes dominant with increasing M as we show below. Compared 
to the inviscid eigenfunction qivs, the viscous eigenfunction q decays due to viscous dissipation 
and thermal diffusion (both of which involve viscosity /i). Under the assumed viscosity-law (O, 
the viscosity of the mean flow increases rapidly with increasing Mach number at all points in the 
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normal direction (see Fig. 1), resulting in a decrease in the viscous eigenfunction q in the same 
limit. To quantify the last statement, we calculate the following for the viscous problem: 

^ El lo Tr qWjt, y)q(^\t, y)dy Ei,k MO)\'e' Jp' \q'i'\y)\'dy 

EJo^Trg(Ot(0,y)g(0(0,y)dy E.MOW lo Wi^vWdy ' 

which is a measure of the collective evolution (growth/decay) of all eigenfunctions. Here, the 
summation index / runs over the perturbation velocity, density and temperature, and the index k 
runs over the selected eigenmodes (see eqn. 19); Im stands for the imaginary part, and Tr stands 
for the trace of the matrix. The initial condition, {k(0)}, is chosen as {1}. It should be pointed 
out that the measure for the collective growth/decay of eigenfunctions, via (1431) . is equivalent to 
probing the energy-norm with an unit weight matrix Ai = I. Our definition simply masks out the 
Mach-number dependence of the energy norm due to the base-state variables in 7V1 (as in the Mack 
energy norm). 

Figure [17] shows the variation of with time for two Mach numbers; other parameters are 
a = and /3 = 1. It is clear from this figure that the (collective) decay-rate of eigenfunctions at 
any time is higher for higher Mach number. Therefore, the viscous eigenfunctions q are small at 
a time t = t^pt or tmax and at a Mach number Mq in comparison with that at any lower value of 
the Mach number M < Mq, leading to a decrease in the energy transfer from the mean flow to 
perturbations with increasing Mach number. This decay can only come from the dependence of 
the viscous eigenfunctions on the shear viscosity (which increases with increasing M). In contrast, 
the inviscid eigenfunctions grow with M as we have pointed out earlier. Therefore, the dijference 
in the variation of Si with M for the inviscid and viscous cases stems from the different variations 
of the respective eigenfunctions with M. 

The reason for the viscous growth curves (in Fig. 16) not coinciding with the inviscid curve 
during the most part of the transient growth can be related to the fact that the viscous eigenfunctions 
undergo sharp decays with time for the present mean flow. (Note that the non-dimensional viscosity 
/io takes a value that is always greater than one in the entire domain, see Fig. Id.) Here, the 
viscous-decay of eigenfunctions with time (as in Fig. 17) competes with the inviscid growth of 
eigenfunctions. Clearly, the decay of Si with time is not negligible compared to the algebraic 
growth even duing the initial times, and this decrease becomes more and more important at later 
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Figure 16: Inviscid energy growth (Giys{t)) and its comparison with viscous growth curves at 
various Reynolds number for M = 2, a = and (3 = 1. The inset shows the variations of Giys{t) 
for three different Mach numbers. 
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Figure 17: Variation of the collective growth/decay of all eigenfunctions, S^f, with time: a = 0, 
/? = 1 and i^e - 2 X 10^ 
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times. Hence, the G(t)-curves for various Re do not coincides with their inviscid counterpart 

Returning to the compressible boundary layers (13], we had noted that the viscous growth 
curves, G{t), coincide with the inviscid curve Givsl^) till they achieve their maxima Gmax- A 
plausible explanation for this could be that the temperature, and hence the viscosity, for boundary 
layers remains constant for more than 99% of the domain (since the variations of mean-fields are 
concentrated within a few displacement thickness while the mean flow extends upto a few hundred 
times of the displacement thickness). Therefore, the viscous-decay of q, compared to the inviscid 
(algebraic) growth, may not be as strong in boundary layers as in the present mean flow. This issue 
needs further investigation. 



6 Summary and Conclusion 

We have investigated the nonmodal transient growth characteristics and the related patterns in 
compressible plane Couette flow of a perfect gas with temperature-dependent viscosity. The mean 
flow consists of a non-uniform shear-field and varying temperature and viscosity along the wall- 
normal direction. For the transient growth analysis, the disturbance size was measured in terms of 
the Mack energy norm [l2l|3l for which the pressure-related energy transfers are zero. The results 
were presented for ranges of Mach number (Af ), Reynolds number (Re) and wavenumbers {a and 
P) for which the flow is asymptotically stable. 

The maximum transient energy growth, Gmax^ is found to increase with increasing Reynolds 
number, as in many incompressible flows, but decreases with increasing Mach number. The opti- 
mal energy growth, Gopt, (i.e. the global maximum of Gmax in the (a, /3)-plane for given Re and 
M) decreases with increasing M. This result is in contrast to that for compressible boundary lay- 
ers [[T2l for which Gopt increases with increasing M. The optimal streamwise wavenumber, a^pt, 
is close to zero (but finite) at M 0, increases with increasing M and reaches a maximum value 
at some value of M, and decreases thereafter. Unlike in incompressible Couette flow [I51|9l, Oopt 
becomes zero at large enough value of M. The optimal spanwise wavenumber, P^pt, also varies 
non-monotonically with M: /5opt decreases first and then increases with increasing M. Optimal 
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velocity patterns (at t = tmax) correspond to pure streamwise vortices for large M, but the mod- 
ulated streamwise vortices are optimal patterns for low-to-moderate values of M. Our result on 
optimal patterns at very high Mach number should be contrasted with that for the incompressible 
Couette flow for which the oblique modes constitute the optimal patterns. 

For the streamwise independent disturbances (a = 0), we have found that the transient energy 
growth does not follow the well-known scaling laws, Gmax ~ R^^ and tmax ~ -Re, of incompress- 
ible shear flows [fT0l l9ll. In contrast, however, these scaling laws are known to hold for compressible 
boundary layers lfT2l . We showed that the invalidity of these scaling laws for the present flow con- 
figuration is tied to the 'dominance' of some terms (related to density and temperature fluctuations 
in the y and z-momentum equations) in the linear stability operator. More specifically, we found 
that the well-known Mack transformation (eqn. [26l) does not make the streamwise-independent 
stability equations independent of the Reynolds number because of the above mentioned dominant 
terms. 

An evolution equation for the perturbation energy has been derived, and various constituent 
energies, that are transferred to perturbations through different physical processes, have been iden- 
tified. We have carried out a detailed nonmodal energy analysis for initial perturbations that yield 
maximum energy growth at a later time. Based on this energy budget analysis, we found that the 
transient energy growth occurs due to the transfer of energy from the mean flow to perturbations 
via an inviscid algebraic instability. We further showed that the decrease of transient growth with 
increasing Mach number is tied to the decrease in the energy transferred from the mean flow {8i) 
in the same limit. Lastly, considering the inviscid limit of stability equations, we found that the 
viscous growth curves follow the inviscid growth curve {Gi^s) only for a very short time. This is 
due to the strong dependence of viscosity on Mach number for the present mean flow, resulting in 
sharp decays of the viscous eigenfunctions with increasing Mach number which is responsible for 
the decrease of Ei in the same limit. 
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